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Distributed  Beamforming  in  Wireless  Relay 
Networks  with  Quantized  Feedback 

Erdem  Koyuncu,  Yindi  Jing,  and  Hamid  Jafarkhani 

Abstract 

This  paper  is  on  quantized  beamforming  in  wireless  amplify-and-forward  (AF)  relay  networks.  We 
use  the  Generalized  Lloyd  Algorithm  (GLA)  to  design  the  quantizer  of  the  feedback  information  and 
specifically  to  optimize  the  bit  error  rate  (BER)  performance  of  the  system.  Achievable  bounds  for 
different  performance  measures  are  derived.  First,  we  analytically  show  that  a  simple  feedback  scheme 
based  on  relay  selection  can  achieve  full  diversity.  Unlike  the  previous  diversity  analysis  on  the  relay 
selection  scheme,  our  analysis  is  not  aided  by  any  approximations  or  modified  forwarding  schemes. 

Then,  for  high-rate  feedback,  we  find  an  upper  bound  on  the  average  signal-to-noise  ratio  (SNR)  loss 
and  show  that  it  decays  at  least  exponentially  with  the  number  of  feedback  bits,  B.  Using  this  result,  we 
also  demonstrate  that  the  capacity  loss  also  decays  at  least  exponentially  with  B.  In  addition,  we  provide 
approximate  upper  and  lower  bounds  on  the  BER,  which  can  be  calculated  numerically.  Simulations  are 
also  provided,  which  confirm  our  analytical  results.  We  observe  that,  for  R  relays,  our  designs  achieve 
full  diversity  when  B  >  log  R  and  a  few  extra  feedback  bits  are  sufficient  for  a  satisfactory  performance 
in  terms  of  the  array  gain.  Simulations  also  show  that  our  approximate  BER  is  a  reliable  estimation  on 
the  actual  BER  for  even  moderate  values  of  B. 

Index  Terms 

Wireless  relay  networks,  beamforming,  vector  quantization,  bit  error  probability,  diversity  order. 

I.  Introduction 

Multiple  antennas  can  improve  the  performance  of  communication  systems  by  introducing 
diversity.  The  diversity  benefit  is  a  result  of  having  statistically  independent  communication 

This  work  was  supported  in  part  by  an  ARO  Multi-University  Research  Initiative  (MURI),  grant  #  W91  lNF-04-1-0224. 

The  authors  are  with  the  Department  of  Electrical  Engineering  and  Computer  Science  at  the  University  of  California,  Irvine, 
CA,  92697-2625.  E-mail:  {ekoyuncu,  yjing,  hamidj}@uci.edu 
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paths  between  the  transmitter  and  the  receiver.  However,  the  use  of  multiple  antennas  may  not 
be  practical  due  to  the  power  and  cost  limitations  of  the  mobile  devices.  Recent  developments 
showed  that  such  difficulties  can  be  overcome  by  the  concept  of  cooperation,  in  which  several 
relays  assist  the  communication  between  the  terminals  of  the  network.  In  this  case,  each  relay 
provides  an  independent  communication  path,  and  this  results  in  a  cooperative  diversity  gain 
[1H3]. 

Many  cooperative  schemes  have  been  proposed  [4]-[7],  among  which  non-regenerative  schemes 
such  as  amplify-and-forward  (AF)  have  especially  become  popular  due  to  their  low  complexity. 
When  channel  state  information  (CSI)  is  not  available  at  the  relay  terminals,  distributed  space- 
time  coding  schemes  have  been  established  [5],  [8],  [9],  which  achieve  full  diversity.  On  the 
other  hand,  with  available  CSI,  each  relay  should  adjust  its  transmit  power  and  its  transmit 
signal  phase  for  a  better  performance.  This  adaptive  network  design  concept  is  analogous  to 
beamforming  in  multiple  antenna  systems,  and  is  thus  called  network  beamforming  [10]. 

The  beamforming  problem  for  networks  with  parallel  relays  is  studied  in  [11]-[14]  for  a  sum- 
power  constraint  on  relays,  and  in  [10],  [12],  for  networks  with  individual  power  constraints. 
Both  works  require  that  each  relay  knows  its  own  channels  perfectly,  in  addition  to  one  or  two 
real  numbers  that  can  be  calculated  by  the  receiver  and  common  for  all  relay  nodes.  A  more 
practical  assumption  is  that  there  is  only  partial  channel  information  at  the  relay  nodes. 

There  are  two  widely  used  types  of  partial  channel  information:  channel  statistics  (means  and 
covariances)  and  quantized  instantaneous  CSI.  The  first  case  was  addressed  in  [15].  In  this  work, 
we  focus  on  the  latter.  In  [16],  outage  minimization  of  single-relay  networks  with  quantized  CSI 
is  performed.  It  is  assumed  that  there  is  a  long-term  power  constraint  on  the  total  power  of  the 
transmitter  and  the  relay.  A  numerical  treatment  for  general  parallel  AF  networks  with  quantized 
CSI  can  be  found  in  [17].  Note  that,  a  special  case  of  quantized  feedback  is  relay  selection  [18]- 
[24],  which  uses  a  fixed  number  of  logi?  feedback  bits  for  a  network  with  R  relays. 

Although  the  concept  of  network  beamforming  with  quantized  channel  information  is  highly 
similar  to  quantized  beamforming  in  multiple-antenna  systems  [25]-[29],  we  cannot  directly 
apply  the  results  on  the  latter.  This  is  due  to  the  existence  of  the  extra  channels  between  the 
transmitter  and  the  relays  that  results  in  noise  amplification  at  the  relay  terminals. 
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In  this  work,  we  analyze  networks  with  any  number  of  relays,  in  which  each  relay  has  its 
own  short-term  power  constraint.  We  assume  full  channel  information  at  the  receiver,  but  each 
relay  only  knows  its  own  receiving  channel  and  a  partial  CSI  provided  by  B  bits  of  feedback. 
The  feedback  bits  are  conveyed  from  the  receiver  to  the  relays  and  they  represent  a  quantized 
network  beamforming  vector.  This  beamforming  vector  is  selected  from  a  codebook  which  is 
known  at  the  relays  and  the  receiver.  Our  performance  measure  is  the  bit  error  rate  (BER). 

The  paper  is  organized  as  follows:  In  Section  II,  we  introduce  our  network  model  and  the 
problem  definition.  In  Section  III,  we  describe  our  codebook  design  strategy  based  on  the 
Generalized  Lloyd  Algorithm  (GLA).  In  Section  IV,  we  investigate  the  achievable  diversity 
of  the  optimal  relay  selection  scheme,  and  show  that  B  feedback  bits  can  achieve  diversity  2^ 
for  B  =  0,  log  2, . . . ,  log  /?.  In  Section  V,  we  provide  an  upper  bound  on  the  average  signal-to- 
noise  ratio  (SNR)  loss  for  high-rate  feedback  and  also  extend  this  result  for  the  capacity  loss.  In 
Section  VI,  we  derive  approximate  upper  and  lower  bounds  on  the  BER,  which  give  an  accurate 
description  of  the  actual  BER  for  moderate  and  large  values  of  B.  The  numerical  results  are 
provided  in  Section  VII.  Einally,  in  Section  VIII,  we  draw  our  major  conclusions.  Some  technical 
proofs  are  provided  in  the  appendices. 

Notation:  ||  ■  ||  indicates  the  2-norm,  (-I-)  is  the  inner  product.  Z"*",  M,  and  C  represent 
the  sets  of  positive  integers,  real  numbers  and  complex  numbers,  respectively.  Eor  2;  G  C, 
|;2|  indicates  the  absolute  value,  2:*  denotes  the  complex  conjugate,  and  y{t{z)  is  the  real  part 
of  ;2.  arg2:  is  the  complex  argument  of  2:  with  0  <  argz  <  2tt.  E[-]  is  the  expected  value  and 
Pr  represents  the  probability.  fx{-)  is  the  probability  density  function  (PDE)  of  the  random 
variable  X.  X  ~  Exp (^2;)  means  that  X  is  an  exponential  random  variable  distributed  as 
fxix)  =  CxGwi~Cxx)  for  a;  >  0  and  fxix)  =  0  for  a;  <  0.  Eor  any  sets  A  and  B,  A  —  B  is 
the  set  of  elements  in  A,  but  not  in  B.  A  c  B  means  ^  is  a  subset  of  B,  and  ^  fi  is  the 
intersection  of  A  and  B.  \A\  is  the  cardinality  of  A  and  A^  =  {(oi, . . . ,  a^)|ai, . . . ,  G  A}, 
r  G  Z+,  is  the  cartesian  power.  Einally,  0  is  the  empty  set,  Q(-)  represents  the  Gaussian  tail 
function,  P(-)  is  the  gamma  function  and  K^{-)  is  the  modified  Bessel  function  of  the  second 
kind  of  order  u. 


II.  Network  Model  and  Problem  Statement 
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A.  System  Model 

The  block  diagram  of  the  system  is  shown  in  Fig.  1.  We  have  a  relay  network  with  one 
transmitter-receiver  pair  and  R  parallel  relays.  We  assume  that  there  is  no  direct  link  between 
the  transmitter  and  the  receiver.  Denote  the  channel  from  the  transmitter  to  the  ith  relay  by  /* 
and  the  channel  from  the  iih  relay  to  the  receiver  by  gi.  Let  h  =  (/i, . . . ,  fn,  gi,  ■  ■  ■  ■,  gn)  be  the 
corresponding  channel  vector.  The  entries  of  h  are  assumed  to  be  independent,  and  each  entry 
is  modeled  as  CJ\f{0, 1),  a  complex  Gaussian  random  variable  with  variance  1/2  per  complex 
dimension. 

Only  the  short-term  power  constraint  is  considered,  which  means  that  there  is  an  upper  bound 
on  the  power  used  at  each  node,  and  for  every  symbol  transmission,  the  power  levels  used  at 
the  transmitter  and  the  ith  relay  are  no  larger  than  Pq  and  Pi,  respectively. 

We  assume  a  quasi-static  channel  model,  in  which  the  channel  realizations  vary  independently 
from  one  frame  to  another,  while  within  each  frame  the  channels  remain  constant.  We  assume 
that  the  receiver  knows  h  and  each  relay  has  perfect  knowledge  of  its  own  receiving  channel, 
namely  the  ith  relay  knows  /*.  Each  relay  also  has  B  bits  of  partial  CSI,  provided  by  the 
receiver  feedback.  The  feedback  channels  are  assumed  to  be  error-free  and  delay-free.  Note  that 
within  each  frame,  the  information  each  relay  receives  from  the  feedback  is  the  same.  With  this 
property,  there  is  a  common  codebook  C  =  |x^, . . .  ,x^|  of  cardinality  M  =  2^,  which  is 
designed  offline  and  known  by  all  the  relay  terminals  and  the  receiver. 

The  feedback  transmission  scheme  operates  as  follows:  For  each  frame,  the  channel  realization 
h  is  quantized  using  a  quantizer  Q  defined  by  the  encoder  and  the  decoder  mappings  S  :  X 

and  V  :  X  ^  C,  where  X  =  {1, . . . ,  M}  represents  the  index  set  for  the  codebook  elements. 
The  encoding  operation  is  performed  at  the  receiver  and  the  feedback  bits  represent  the  encoder 
output.  Then,  each  relay  decodes  the  feedback  bits  using  V,  to  find  the  corresponding  codebook 
element.  If  i^(h)  =  m,  we  have  V{m)  =  x"^  G  C.  For  the  rest  of  this  paper,  we  will  use  the 
well-known  representation  Q(h)  =  V{S{\i))  as  the  combined  effect  of  the  mappings  S  and  V. 
Therefore,  Q  ;  C  and  Q(h)  =  x™,  for  some  m  G  X  and  x™  G  C. 
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B.  Transmission  Scheme  and  Performance  Measure 

We  consider  a  two-step  AF  protocol  [10].  During  the  first  step,  the  transmitter  selects  a  symbol 
s  from  a  constellation  S,  and  sends  s/^s.  We  normalize  s  as  E[|sp]  =  1.  Therefore,  the  average 
power  used  at  the  transmitter  is  Pq.  During  this  step,  there  is  no  reception  at  the  receiver,  but  the 
Ah  relay  receives  =  s/^fiS  -f  z/*.  The  noises  z/*  are  independent  and  modeled  as  CA/’(0, 1). 

One  of  the  key  issues  in  cooperative  network  design  is  the  relay  operation.  In  this  work,  we  use 
AF  relays.  However,  unlike  [4]-[6],  we  incorporate  the  idea  of  network  beamforming  [10],  by 
allowing  each  relay  to  adaptively  adjust  its  transmit  power  and  phase  according  to  the  feedback 
information.  Thus,  each  element  of  the  codebook  C  corresponds  to  a  beamforming  vector. 

Assume  that  the  receiver  feeds  back  the  index  m.  Then,  the  relays  use  the  beamforming  vector 
=  {xf, . . .  ,xff).  In  the  second  step,  the  Ah  relay  sends  =  xf^/plrt  where  pi  = 

The  average  transmit  power  of  the  Ah  relay  can  be  calculated  as  \xf\‘^Pi,  where  <  \xf\  <  1 
is  a  result  of  the  short  term  power  constraint. 

After  these  two  steps  of  transmission,  the  received  signal  at  the  receiver  can  be  expressed  as: 

y  =  fe  S  +  J2 

\j=l  /  i=l 


where  oo  is  the  noise  at  the  receiver,  which  is  also  distributed  as  CJ\f{0, 1)  and  independent  from 
z/j.  The  received  SNR  can  be  given  as: 


7(h;  Q,C) 


Po 


l  +  Ell\xT\^\9^?P^ 


(1) 


for  channel  state  h  with  the  quantizer  Q  and  the  beamforming  codebook  C.  We  will  also 
frequently  use  the  notation  7(x,  h)  to  describe  the  received  SNR  for  a  specific  beamforming 
vector  X.  Note  that  given  C,  7(h;  Q,  C)  =  7(x,  h)  whenever  x  =  Q(h). 

In  this  paper,  our  performance  measure  is  the  BER,  which  can  be  simply  expressed  as  an 
expectation  over  the  channel  state  vector: 


where  (3  is  a  constant  depending  on  the  constellation  S. 
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C.  Problem  Statement 

Our  network  beamforming  design  problem  is  the  joint  optimization  of  the  quantizer  and  the 
quantization  eodebook,  given  a  fixed  number  of  feedbaek  bits  B,  sueh  that  the  BER  is  minimized. 
The  quantizer  design  is  straightforward: 

Lemma  1:  For  any  beamforming  eodebook  C  =  |x^, . . . ,  x^|,  the  optimal  quantizer  Q*  to 
minimize  the  BER  is  given  by: 

Q*(h)  =  argmax7(x,  h).  (3) 

xSC 

Proof:  Consider  a  fixed  eodebook  C  and  an  arbitrary  Q!  :  C.  For  any  h,  we  have 

7(h;  Q*,C)  >  7(h;  Q',C)  and  the  result  follows  from  the  monotonieity  of  the  Q  function.  ■ 
In  other  words,  given  C,  the  optimal  quantizer  chooses  the  SNR  maximizing  beamforming 
vector  in  C.  From  now  on,  we  fix  the  quantizer  operation  to  be  Q*.  Therefore,  the  eodebook 
uniquely  determines  the  performance  of  our  feedback  system.  Before  introducing  the  eodebook 
design  methods,  we  give  a  necessary  condition  for  a  eodebook  to  be  optimal. 

When  there  is  no  limit  on  the  number  of  feedback  bits,  the  relays  can  always  use  the  optimal 
beamforming  vector,  which  was  derived  in  [10].  In  this  case,  regardless  of  the  channel  state,  at 
least  one  relay  should  transmit  with  full  power.  This  requirement  can  be  represented  by  the  set: 

A"  =  {x  :  X  G  and  ||x||oo  =  1},  (4) 

where  ||x||oo  —  max*  \xi\.  For  our  network  model  with  finite  feedback  bits,  we  show  that  each 
vector  in  an  optimal  eodebook  should  also  be  an  element  of  the  set  X\ 

Lemma  2:  If  C*  is  optimal,  then  C*  (Z  X. 

Proof:  Suppose  an  optimal  eodebook  C  with  C  —  X  ^  exists.  So,  there  exists  x^  G  C 
where  \x\\  <  1,  Vi  G  {!,. . .  ,i?}.  Define  7  =  max*  \x\\  and  let  y  =  ^x^.  Replace  yf  in  C 
with  y,  which  results  in  a  new  eodebook  C  =  {x\  . . . ,  x^“^,  y,  x^+^, . . . ,  x^}.  From  (1),  it  is 
easy  to  see  that  7(y,h)  >  7(x^,h),  Vh  G  Therefore,  'P{E\Q*,C)  <  'P{E\Q*,C')  and  C  is 
suboptimal.  This  contradicts  the  optimality  of  C'.  ■ 

This  lemma  simply  states  that  at  least  one  element  of  every  beamforming  vector  in  an  optimal 
eodebook  should  have  unit  norm.  Therefore,  the  set  X  completely  defines  the  feasible  set  for  the 


7 


codebook  vectors.  Note  that,  for  networks  with  a  sum-power  constraint  on  relays,  the  feasible 
set  of  codebook  vectors  can  similarly  be  defined  by  the  set  {x  ;  x  G  and  ||x|p  =  1}. 


III.  CoDEBOOK  Design:  The  Generalized  Lloyd  Algorithm 

We  use  the  Generalized  Lloyd  Algorithm  (GLA)  to  design  optimal  codebooks  of  an  arbitrary 
size.  The  algorithm  is  based  on  two  optimality  conditions,  which  can  be  re-stated  as  follows: 

For  a  given  codebook  {x"^  :  m  =  1, . . . ,  M},  the  optimal  partition  of  the  channel  state  space 
is  given  by: 

=  {h  :  7(x”",  h)  >  7(x’",h),Vn  7^  m},  (5) 


and  for  a  given  partition  {IZm  :  m  =  1, . . . ,  M},  the  optimal  codebook  vectors  satisfy: 


=  argminEh 

xGA:’ 


Q  (\/2/37(x,h)]  I  h  e 


(6) 


In  general,  the  algorithm  is  initialized  with  a  random  choice  of  codebook  vectors.  Then,  the 
optimality  conditions  are  iterated  until  convergence  is  achieved.  Note  that  the  GLA  does  not 
guarantee  global  optimality.  However,  at  each  iteration,  the  algorithm  provides  an  improved 
codebook  design  and  each  new  codebook  results  in  a  lower  BER.  Therefore,  it  is  guaranteed 
that  the  algorithm  converges  and  the  resulting  codebook  achieves  at  least  the  BER  performance 
of  the  initialization  codebook. 

The  GLA  can  also  be  used  for  designing  beamforming  codebooks  for  networks  with  a  sum- 
power  constraint  on  the  relays.  The  only  difference  is  in  the  shape  of  the  feasible  set.  Moreover, 
all  the  diversity  and  scaling  results  we  will  obtain  for  networks  with  individual  power  constraints 
apply  to  networks  with  a  sum-power  constraint  as  well. 


IV.  Diversity  Analysis  oe  Relay  Selection  Scheme 

The  maximal  spatial  diversity  of  networks  with  R  relays  is  R.  In  this  section,  we  will  prove 
that  a  simple  feedback  scheme  based  on  relay  selection  achieves  this  full  diversity.  Note  that  we 
can  also  improve  the  relay  selection  codebook  by  GLA.  This  way,  we  make  sure  that  the  designed 
codebook  provides  full  diversity  and  also  better  performance  compared  to  relay  selection. 

For  the  relay  selection  scheme,  only  one  of  the  R  relays  is  allowed  to  cooperate  for  a  given 


constant  fading  block.  The  corresponding  codebook  is  given  by  =  {e™  :  m  =  1, . . . ,  M}, 
where  e™  =  1  if  i  =  m,  and  e™  =  0  for  i  7^  m.  Using  (1)  and  (3),  the  received  SNR  is 


7(h;  Q*,Ce)  =  ,  max 


PoPi  \fi9i 


(7) 


A  BER  analysis  for  this  scheme  was  previously  performed  in  [20]-[22],  which  however,  use 
approximations  on  (7).  In  this  paper,  we  find  a  rigorous  upper  bound  on  the  BER. 

Theorem  1:  Let  Pi  =  XiPq,  where  >  0  are  constants.  Define  A  =  maxj{l  +  which 

is  also  a  constant.  Then,  the  BER  of  the  relay  selection  scheme  is  bounded  by: 

PT{E\Q*,Ce)  <  RPq^X^^  E  ~  (8) 

Proof:  See  Appendix  I.  ■ 

Notice  that  the  exponent  of  Pq  in  (8)  is  —R.  Thus,  the  relay  selection  scheme  achieves  the  full 
diversity  of  R  with  M  =  R.  Similarly,  diversity  order  M  is  achievable  for  M  <  R,  simply  by 
considering  the  selection  scheme  for  any  fixed  M  of  the  R  relays  and  disregarding  the  others. 

As  we  have  discussed  before,  we  can  design  a  new  quantization  codebook  to  achieve  better 
array  gain  compared  to  the  relay  selection  codebook  using  the  GLA.  Since  the  new  codebook 
also  provides  the  same  full  diversity  order  property  as  the  relay  selection  scheme,  its  superior 
average  SNR  results  in  a  better  array  gain. 


V.  Average  SNR  Loss  with  Quantized  Eeedback 

In  this  section,  the  high-rate  feedback  scenario  is  considered.  We  introduce  a  systematic 
codebook  design  Cu,  which  also  incorporates  the  relay  selection  vectors.  Thus,  the  resulting 
design  provides  full  diversity  order.  In  addition,  we  show  that  the  average  SNR  loss  of  the 
design  decays  exponentially  with  the  number  of  feedback  bits,  which  indicates  an  increasing 
array  gain.  We  also  use  this  result  to  find  an  upper  bound  on  the  capacity  loss. 


A.  Codebook  Design 

A  beamforming  vector  contains  both  power  and  phase  information.  Eor  the  design  of  C„,  we 
treat  them  separately  and  construct  two  different  codebooks  to  represent  each.  Then,  the  two 
codebooks  are  combined  to  construct  Cu- 


Quantized  phase  codebook:  We  write  the  reeeived  SNR  as  7(x,  h)  =  Po 


hiPIffiPPi 


where  Ai  =  0  and  Aj  =  argXj  —  argxi,  i  =  2, . . . ,  R.  Without  loss  of  generality,  we  only  need 
to  quantize  A2, . . . ,  A/j.  For  eaeh  A*  e  [0,  2tt),  we  use  a  scalar  quantizer  defined  hy  the  centroids 
U  =  |o,  271^, . . . ,  27r^^|,  L  >  2,  L  e  Then,  we  use  the  set  of  vectors  Cphase  —  to 
quantize  the  phase. 

Quantized  power  codebook:  Similarly,  we  define  V  =  |o,  ^  ij,  K  >  2,  K  e 

and  consider  the  set  V^.  As  a  result  of  Lemma  2,  we  remove  the  vectors  in  hut  not  in 
X,  which  results  in  Cpower  —  fl  A  =  —  {V  —  {1})^.  Fig.  2  depicts  an  example  for  the 

construction  of  the  quantized  power  codebook. 

Therefore,  we  have  |C„|  =  | Cphase || Cpower |  =  \k^  -  {K  - 


B.  Average  SNR  Loss 

The  pointwise  SNR  loss  of  C„  can  be  defined  as  P(h)  A  ^(h;  Q*,  X)  —  7(h;  Q*,  Cu),  which 
is  the  difference  between  the  received  SNRs  with  perfect  feedback  and  with  Cu-^ 

Theorem  2:  The  average  SNR  loss  with  quantized  feedback  decays  at  least  linearly  as  the 
number  of  quantization  vectors  M  =  2^,  or  exponentially  as  the  number  of  feedback  bits  B: 

Eh[P(h)]  </i2“^(^,  (9) 

where  p  is  independent  of  B. 

Proof:  See  Appendix  II.  ■ 

Using  Theorem  2,  we  can  find  a  similar  result  for  capacity  loss  with  quantized  feedback. 


C.  Capacity  Loss 

Let  C  and  Cu  denote  the  ergodic  capacities  with  perfect  feedback  and  with  the  codebook 
Cu,  respectively.  Note  that  U  =  Eh  [log(l  +  7(h;  Q* ,  X))]  and  Cu  =  Eh  [log(l  +  7(h;  Q*,  Cu))], 
both  in  nats  per  channel  use.  Then,  the  capacity  loss  Ci,  =  C  —  Cu  can  be  bounded  as: 


Cr  =  Ev 


log  1  + 


D{h) 


1  +  7(h;  Q*  Cu 


<  Eh 


D{h) 


1  +  7(h;  Q*  Cu) 


<Eh[D{h)]  <  p2  (10) 


'By  Q*,  we  mean  the  optimal  quantizer  with  a  particular  codebook,  thus  the  two  Q*s  in  the  SNR  loss  are  actually  different. 
However,  with  an  abuse  of  notation,  we  omit  the  dependence  of  Q*  to  the  codebook,  as  it  is  clear  from  the  context. 
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where  the  first  inequality  is  a  result  of  the  fact  that  log(l  +  x)  <  a;  for  a;  >  0. 

We  note  that  the  constant  /a  will  not  result  in  a  tight  upper  bound  for  moderate  M,  therefore 
(9)  and  (10)  are  useful  especially  with  high  rate  feedback.  However,  these  results  show  how  the 
average  SNR  loss  and  the  capacity  loss  scale  with  the  number  of  feedback  bits. 

VI.  Bit  Error  Probability  With  Quantized  Feedback 

In  this  section,  we  analyze  the  BER  of  the  network.  One  big  challenge  for  the  BER  analysis 
comes  from  the  shape  of  the  feasible  set  of  the  codebook  vectors  X.  Even  in  the  case  of 
perfect  channel  knowledge  at  the  relay  terminals,  this  leads  to  complicated  expressions  for  the 
optimal  beamforming  vector,  and  consequently,  great  difficulties  in  the  derivation  of  the  system 
performance.  To  overcome  these  problems,  we  bound  the  BER  with  feasible  set  X  by  BERs 
with  feasible  sets  that  are  much  easier  to  deal  with. 

A.  Performance  Bounds  Using  Sum-Power  Constraints 

For  any  beamfoming  vector  x,  first,  we  consider  the  restriction  ||xp  =  1,  which  just  implies 
a  sum-power  constraint  on  relays.  The  corresponding  feasible  set  for  the  codebook  vectors  is 
A"™  =  {x  ;  X  G  and  ||x|p  =  1}.  Similary,  we  define  =  {x  ;  x  G  and  ||x|p  =  R}. 
Note  that  for  any  codebook  C  A“  (or  C  the  optimal  encoder  is  also  Q*. 

Lemma  3:  Let  C™*  C  A™,  C*  C  A  and  C  A°“*  be  the  optimal  codebooks  of  cardinality 
M  for  their  respective  feasible  sets.  Then,  Pr(i?|Q*,C°“**)  <  Pi{E\Q* ^C*)  <  Pr(i?| Q*, C“*). 

Proof:  First,  we  prove  Pi{E\Q* ,C*)  <  Pr(i?| Q*, C“*).  For  any  G  choose  r]m  = 
maxj  \yf\,  and  construct  a  new  codebook  C  =  {^y^,  m  =  1, . . . ,  M}.  Note  that  C  G  A  and 
^  7(y™')h),  for  any  h  G  and  m  G  M}.  Therefore,  Pi{E\Q*  ,C*)  < 

Pr(i?|Q*,C)  <  Pr(i?| Q*, C™*),  as  C*  is  optimal  for  the  region  A. 

Similary,  using  C* ,  we  can  construct  =  {||^^x™  :  ^  C* ^  m  =  Since 

||x’”||  <  i?  as  x”^  G  C*  C  A,  7(||^^x’”,  h)  >  7(x™',  h)  for  any  h  G  and  m  G  {1, . . . ,  M}. 
Therefore,  Pt{E\Q* <  Pr(E|Q*,C°“*)  <  Pt{E\Q*,C*),  as  is  optimal  for  A°"P  ■ 

Therefore,  the  achievable  performance  with  the  feasible  set  A  is  bounded  by  the  achievable 
performance  with  the  feasible  sets  A“  and  A°“\  for  any  number  of  feedback  bits. 
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B.  BER  Analysis  for  Sum-Power  Constrained  Relays  with  Quantized  Feedback 


We  first  consider  the  BER  with  C  For  simplicity,  we  assume  Pi  =  ...  =  Pr  =  P. 
For  a  particular  channel  state  h,  let  y"^  =  Q*(h)  G  C™  for  some  m  G  {1, . . . ,  M}.  Note  that, 
h  G  TZrn  by  the  optimality  conditions  stated  in  (3)  and  (5).  In  other  words,  the  quantization  of  the 
beamforming  vector  also  implies  a  partition  of  the  channel  state  space  into  M  quantization 
cells.  We  consider  the  distribution  of  the  received  SNR  given  h  G  TZm- 

Fet  bi  =  and  b  =  (6i, . . . ,  hn).  We  re-write  the  received  SNR  given  by  (1)  as: 

,  m  _  Po|(b,y"-)p  _Po||bf|('b,^)P_  ^ 

'  1  I  V-R  \yT?\9i?P  i,srR\yTfMZ.  ^o|\u,y  /|  2^  \yf?\gj?P  ^  ’ 

■^^*=1  l  +  l/ipPo  l+|/i|2Po  ^  +  2^j  =  l  l  +  |/^.|2po 


l+l/iP^’o 


where  for  any  complex  vector  z,  we  define  "z"  =  z/||z||.  The  second  inequality  follows  from  the 
fact  that  ||y"^|p  =  1.  To  calculate  the  average  BER,  we  need  the  density  function  of  the  received 
SNR.  However,  the  signal  and  the  noise  powers  of  all  relays  are  coupled  in  the  received  SNR 
formula  given  by  (11),  which  makes  the  density  function  calculation  intractable,  if  not  possible. 
For  analytical  tractability,  we  seek  a  decoupled  approximation  of  (11): 


7(y^h)^PoKb,3^)rE 


\9i?P 

1+I7PP) _ 

'■ — '  1  I  IffipP  Y^P 
=1  -L  +  l^f^\2p^2^j=l  \yj\ 


=  l{b,y”)rr”(h), 


where  T'»(h)  A 

D  l'U™PlO'PP 

the  denominator  l  +  J2j=i  1+7  I 


.  We  can  see  that  for  the  ith  term  in  (11),  we  approximate 


by  1  + 


Id^P  Y^P 


Ef=i  bf  P  =  1+  iTIOT-  In  o*er  words,  for 


the  ith  term,  we  approximate  the  effect  of  noises  from  all  relays  by  the  noise  from  the  ith  relay 
only,  while  preserving  the  noise  power.  Notice  that  the  average  value  of  the  two  are  the  same. 
Thus,  we  have  an  unbiased  approximation  for  the  noise  effect.  In  fact,  T™(h)  is  the  maximal 
received  SNR  for  any  channel  state  h  [11,  Eq.  9],  i.e.,  T“(h)  =  7(h; 

As  T“(h)  is  independent  of  the  quantization  cell  TZm,  the  inner  product  |(b,y™)p  charac¬ 
terizes  the  boundaries  of  the  partitions  and  the  performance  loss  due  to  quantization.  Therefore, 


it  is  sufficient  to  consider  the  density  of  |  ( b  ,  y"^)  p  given  h  G  TZm- 

For  high  SNR,  the  components  of  b  can  be  approximated  as  bi  =  ~  ~ 

\gi\  exp(j(arg/j  -|-  arg^fj)),/^.  Note  that  this  approximation  is  valid  when  |/jp  X.  For  any 
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finite  Pq,  since  |/jp  has  an  exponential  distribution,  there  is  always  a  non-zero  probability  that 
I /ip  3>  is  not  true.  However  as  Pq  increases  to  inbnity,  this  probability  decreases  to  zero. 
Thus,  we  can  say  that  this  approximation  is  valid  almost  surely  for  high  Pq.  Hence,  if  we 
dehne  b'  =  (6'^, . . . ,  6/.),  then  we  have  |  ( b  ,  y”^)  p  ^  |  (b' ,  y”^)  p.  Together  with  (12),  the  actual 
quantization  cells  TZrn  are  approximated  as 


«  1?-'™  =  {h  :  l(t7,?^)P  >  |(l?,r>l",Vn  >n}  . 

In  [26],  the  conditional  density  of  Vt'  =  |(b',y"^)P  given  h  G  TZ'^  has  been  considered  for 
beamforming  in  MISO  systems  with  quantized  feedback.  It  has  been  shown  that  the  random 
variable  has  a  beta  distribution  with  parameters  1  and  R  —  1,  and  for  large  B,  the  conditional 
PDF  of  Vt'  given  h  G  can  be  approximated  by 


/Q/(a;|h  G  7^^) 


0,  otherwise. 


(13) 


Following  the  same  arguments,  we  approximate  the  conditional  PDF  of  D  =  |(b,y'")p  given 
h  G  TZm  by  (13),  for  large  Pq  and  B  as  well. 

Now  let  us  work  on  T“,  which  is  actually  the  sum  of  R  independent  random  variables.  The 
density  function  for  is  simply  the  convolution  of  the  density  functions  that  we  have  derived 
for  the  proof  of  Theorem  1  (See  (18)  in  Appendix  I).  Furthermore,  as  the  norm  and  the  phase  of 
the  components  of  bi  are  independent,  b  is  independent  of  T“.  Thus,  and  D  are  independent. 
Hence,  we  can  combine  (12)  and  (13)  to  find  the  approximated  BER  as: 


=  E. 


Q 

—  En 

EyIii 

Q 

y2/)DT“ 

D  =  a; 

rl  r  roo 

=  2^{R—l)  /  /  Q,[2(3uv]  fyin{v)dv 


fn{uj)duj, 


(14) 


which  can  easily  be  calculated  numerically. 

A  similar  analysis  can  be  carried  out  for  Pr(E';  Then,  it  is  easy  to  show  that 


rpOUt  A  TT' 
H/  — 


Q 


R 


,  where  ^  ^ 


hl^RPoP 


(15) 
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According  to  the  heuristic  discussions  on  the  validity  of  these  approximations,  we  expect  them 
to  work  well  for  large  B  and  high  SNR.  Actually,  the  simulations  in  Section- VII  show  that  they 
are  even  valid  for  moderate  B  and  low  SNR. 

It  is  also  worth  to  point  out  that,  actually  (14)  is  an  accurate  approximation  on  the  BER  for 
networks  with  a  sum-power  constraint  on  the  relays.  In  Section- VII,  we  also  show  the  simulated 
BER  of  networks  with  a  sum-power  constraint  and  compare  this  with  the  approximated  BER. 

VII.  Simulation  Results 

In  this  section,  numerical  results  are  provided.  Unless  otherwise  specified,  we  assume  an  equal 
power  constraint  on  both  the  transmitter  and  the  relays.  In  other  words,  Pq  =  Pi  =  ■  ■  ■  =  Pr  =  P. 
The  horizontal  axis  of  the  figures  indicates  P,  and  the  vertical  axis  represents  the  average  BER. 

Eig.  3  shows  the  simulation  results  for  networks  with  individual  power  constraints  on  relays. 
The  GLA  is  used  to  design  the  quantizer  codebook  for  the  feasible  set  X.  To  avoid  bad  suboptimal 
codebook  designs,  we  also  use  simulated  annealing  [30].  We  can  see  that  the  relay  selection 
achieves  full  diversity  order  in  all  cases.  However,  with  the  same  number  of  feedback  bits,  the 
relay  selection  scheme  is  suboptimal  compared  to  network  beamforming  using  the  GLA.  This 
can  be  observed  for  networks  with  2  and  4  relays,  and  1  and  2  feedback  bits,  respectively. 
Also  notice  that  a  few  feedback  bits  can  result  in  a  very  close  performance  to  the  optimal  case 
(cx)  bits).  At  a  BER  of  10“^,  the  network  with  3  relays  and  4  bits/8  bits  feedback  bits  is  only 
1.5dB/0.5dB  worse  than  the  optimal. 

The  GLA  can  also  be  used  to  design  codebooks  for  networks  with  unequal  individual  power 
constraints.  In  Eig.  4,  we  show  the  performance  results  for  a  2-reIay  network.  The  power 
constraints  of  the  source  and  the  relays  are  chosen  as  Pq  =  Pi  =  P  and  P2  =  f  •  In  Fig. 
5,  we  show  the  performance  of  a  3 -relay  network  in  which  Pq  =  P  and  Pi  =  P2  =  Pq  = 
either  case,  both  the  relay  selection  scheme  and  the  feedback  schemes  achieve  full  diversity. 

In  Eig.  6,  we  show  the  performance  of  networks  with  a  sum-power  constraint  on  relays.  The 
reasons  for  this  are  two-fold.  Eirst,  we  show  that  the  GLA  can  also  be  used  to  design  codebooks 
for  networks  with  a  sum-power  constraint  on  relays,  and  second,  we  show  the  validity  of  our 
BER  approximations  in  (14)  and  (15).  Eirst,  we  can  see  that  network  beamforming  using  the 
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GLA  works  for  this  network  model  as  well.  Full  diversity  is  aehieved  for  2,  4,  and  8  bits  of 
feedbaek  in  3-relay  networks.  Compared  to  the  infinite  feedback  bits  case,  network  beamforming 
with  8  feedback  bits  is  less  than  0.5  dB  worse.  In  Fig.  6,  we  also  compare  the  simulated  BER  of 
the  GLA  designs  with  (14),  for  networks  with  3  relays.  We  can  observe  that  the  approximation 
provides  a  lower  bound  on  the  BER  in  all  of  the  cases.  This  is  intuitively  clear,  as  similar  to 
[26],  the  random  variable  Q  overestimates  the  distribution  of  |  (  b  ,  y"^)  p  within  each  quantization 
cell.  Moreover,  T“(h)  is  the  maximal  received  SNR  for  any  channel  state  h,  and  as  5  — >  cx), 
we  have  |(b,y”^)P  ^  1,  and  the  analysis  remains  consistent  with  the  actual  performance.  The 
approximated  BERs  are  accurate  even  when  B  is  small,  within  1  dB  for  all  of  the  cases. 

Einally,  we  compare  the  bounds  in  (14)  and  (15)  with  the  simulated  BER  for  networks  with 
individual  power  constraints  in  Eig.  7.  In  both  cases,  the  upper  bound  and  the  lower  bound  is 
within  1  dB  and  1.5  dB  of  the  actual  BER,  respectively.  Notice  that  the  performance  of  the 
lower  bound  is  better  than  the  performance  with  oo  bits  feedback.  The  reason  is  that  the  lower 
bound  corresponds  to  the  feasible  set  which  can  provide  a  higher  array  gain  compared 
to  the  actual  feasible  set  X.  On  the  other  hand,  as  the  performance  with  oo  bits  is  not  known 
analytically,  the  upper  and  lower  bound  gives  valuable  information  about  the  performance  of  our 
limited  feedback  scheme.  We  also  note  that,  as  a  result  of  the  nature  of  our  approximations,  the 
upper  bound  may  not  work  well  when  B  is  small. 

VIII.  Conclusion 

We  study  quantized  beamforming  schemes  for  wireless  amplify- and-forward  relay  networks. 
The  Generalized  Lloyd  Algorithm  is  used  for  the  codebook  design  procedure.  We  rigorously 
analyzed  the  achievable  performance  in  terms  of  the  diversity  order,  the  average  SNR  loss  and 
the  capacity  loss.  We  also  provided  approximate  upper  and  lower  bounds  on  the  bit  error  rate, 
which  can  be  calculated  numerically.  Our  analytical  results  were  verified  by  simulations.  It 
is  observed  that  our  design  methods  achieve  full  diversity  order  and  a  few  feedback  bits  is 
sufficient  to  achieve  a  satisfactory  performance  in  terms  of  the  array  gain.  We  also  observed  that 
our  approximate  BER  analysis  provides  an  accurate  characterization  of  the  performance,  even 
for  moderate  values  of  B. 


Appendix  I 
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Proof  of  Theorem  1 


For  simplicity  of  notations,  define  the  random  variables  Xi  =  and  F*  = 

>  0,  Vi,  whieh  implies  Xi  ~  Exp(Cj)  and  Yi  ~  Exp(,^j).  Note  that  (8)  will  follow  from 
the  substitutions  C*  =  1  and  =  A“\  Now,  the  relay  selection  function  can  be  expressed  as 
Hi  =  ^  1, . . . ,  -R,  where  7  =  Pq-  Thus,  is  a  random  variable  whieh  is  also  the 

received  SNR  given  that  the  ith  relay  is  seleeted. 

Lemma  4:  The  eumulative  distribution  function  of  Hi  is  given  by: 


Pr(P,  <  7)  = 


l-exp(^-^{Ci  +  ^i)^  kKi{k)  ,  7>0, 

0,  7  <  0. 


where  k  = 


/4C»67(1  +  ^ 

/y2 


(16) 


Proof:  Let  X  ~  Exp(C)  and  Y  ~  Exp(,^).  Define  H  =  i^^x+'yY’  where  7  is  a  constant. 
Then,  Pr(P  <  7)  =  Pr  (i+ff+^y  <  7)  =  /o“  Pr(7^a;F  -  77F  <  7  +  77a;|X  =  x)fx{x)dx. 
Note  that  Pr(7^a;F  —  77F  <  7  +  77a;|X  =  x,  X  <  7/7)  =  1.  Therefore, 

Pr(P  <  7)  =  fx{x)dx+  r  Pr  (y  <  ^±^|X  =  x]fx{x)dx.  (17) 
Jo  77/7  \  Yx  -  77  / 

The  first  integral  is  obvious  sinee  fx{x)  =  For  the  seeond  integral,  we  substitute  u  = 

a;  —  7/7  and  use  Pr(F  <  i/)  =  1  —  for  1/  >  0.  Thus, 


Pr(P  <  7) 


-C^ 
e  ~i 


+  Ce 


-C5 


1  —  exp 


/  ^7(1  +  7) 
\  u  Y 


e  ^'^du, 


which  can  be  evaluated  using  JY  exp  ^  —  (Jx'^  dx  =  for  a,/)  >  0  [31,  Eq. 

3.324.1].  Then,  the  result  follows  after  some  simplifieations.  The  discontinuity  of  the  function 
Ki{k)  at  n  =  0  for  7  =  0  can  be  avoided  by  a  direet  evaluation  of  (17).  ■ 

It  is  now  easy  to  ealeulate  the  PDF  of  Hi  using  K[{x)  =  —Ko{x)  —x~^Ki{x)  [32,  Eq.  9.6.26] 
and  differentiating  (16)  with  respeet  to  7. 


-7, 


/Hi(7)=exp  — (Ci  +  6)  7  [N7(Ci +  6)^1  («)  +  2C*6^o(n)  (1  +  27)]. 


7 


(18) 


The  reeeived  SNR  for  the  relay  seleetion  seheme  is  given  by  the  random  variable  H  = 
maXi  Hi{Y-  Sinee  ^*(7)  are  independent. 


R 


hii)  =  Y. 


i=l 


R 


n  PrWSi) 


fnAl) 


(19) 
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by  the  rule  for  the  derivative  of  a  product.  Now,  the  BER  Pr(i?)  can  be  expressed  as 


Ft{E)  =  [q  (^w)]  =  Mm  ym)  di.  (20) 

We  want  to  find  an  upper  bound  for  (20).  As  an  intermediate  step,  we  have: 

Lemma  5:  The  following  bounds  apply  for  K^{x),  v  G  z+,  'ix  >  0,  x  G  M: 

M{x)  <  rm{v)x-\  (21) 

xKm  >  e-^  (22) 

Km  <  m.  (23) 

Proof:  We  start  with  (21),  and  use  an  integral  equality  for  the  function  K^{x)  [32, 
Eq.  9.6.25]:  K,(x)  =  cos(xt)X  +  /“(t^  +  1)— Idi  = 

2^~^T{u)x~’^,  where  [31,  Eq.  3.241]  is  used  to  evaluate  the  last  integral. 

Eor  (22),  we  use  Ki{x)  =  x  /“  —  Idt  [32,  Eq.  9.6.23].  As  —  1  >  t  —  1  for  f  >  1, 

KM)  >  e-"*  (t  -l)dt  =  x^-^  = 

Eor  (23),  we  use  the  following  recurrence  relation  [32,  Eq.  9.6.26]:  ~  = 

KCM)i  dlv  —  Substituting  o  =  1,  Kq{x)  =  K2(x)  —  2x~^Ki{x)  <  2x~‘^  —  2x~‘^e~^  < 

2x~^,  where  for  the  first  inequality  we  have  used  (21)  and  (22)  to  bound  K2{x)  and  —Ki{x), 
respectively.  The  last  inequality  follows  from  1  —  <  a;  for  a;  >  0.  ■ 

Now,  we  bound  fuM)  in  (18)  as 


/h,(7)<-  c^+^^+2^/a^yyl= ,  (24) 

^  I.  V7(l  +  7). 

where  we  used  exp  (^(0  +  6))  <  nnd  (21)  and  (23)  to  bound  kKi{k)  and  Ko{k).  Similary, 
for  Pr(i;fj  <  7)  given  in  Lemma  4,  we  have 

1  " 

Fr{Hi  <  7)  <  -  7  (Ci  +  6)  +  2v^\/7(1  +  7)  •  (25) 

Substituting  (24)  and  (25)  to  (19),  we  have  the  representation  fnif)  <  'I' (7)  7”^^  where 
4/  (7)  is  independent  of  7.  Next,  we  show  that  4/  (7)  is  bounded  by  a  polynomial  of  degree 
{R  —  1)  of  7  for  R>2.  Let  a  =  max*  (Ci  +  b  =  max*  Then,  we  have 


/ - r 

'J' (7)  < -R  07  +  6^7(1  +  7)  a  +  b- 

L  -’Ll 


1  +  27 
^7(1  +  7) 
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<R  a\J^{l  +  7)  +  h^^il  +  7)  a  +  h^ 

=  /?  (a  +  b)^~^  7^^  (1  +  7)  ~  (7(0  +  2b)  +  6) 

<  2i?  (a  +  &)^  7^  (1  +  7)  ^  <  2i?  (a  +  6)^  (1  +  7)^“^ .  (26) 


Now  let  A  =  maxj  |v^  +  v^}  =  (a  +  &)2.  Combining  (20)  and  (26), 
Pr(E)  <  7-''  ^  (7)  Q  (y^) 


<  r  (1  +  7)^“^  e-d^d^ 

Jo 

fR-l\  /•“  .  c  1 


=  Y1 


=  r^Rx^^  Y. 


which  concludes  the  proof. 


n=0  \ 

^  (R-l 


n  Jo 


7^-i-”e-^^d7 


V{R  -  =  r^RX^^  Y  ^ 

n=0  \  ^ 


Appendix  II 


Proof  of  Theorem  2 


First,  we  will  introduce  a  suboptimal  quantizer  C„  to  simplify  the  proof.  We 

similarly  define  the  pointwise  SNR  loss  for  C„  with  the  quantizer  Qu  as  iA„(h)  =  7(h;  q*,a:)- 
7(h;  Qu,Cu).  Now,  since  7(h;  Qu,Cu)  <  7(h;  Q*,C„),  we  also  have  iA(h)  <  Du{h),  Vh  e 
which  implies  Eh  [C*(h)]  <  Eh  [iA„(h)].  Therefore,  it  suffices  to  consider  the  average  SNR  loss 
for  the  system  with  codebook  Cu  and  the  quantizer  Qu. 

Let  y  denote  the  optimal  beamforming  vector  for  some  channel  state  h.  By  definition,  we 
have  |r/g|  =  1  for  some  q  G  {1, . . . ,  i?}.  Let  argr/i  =  0,  without  loss  of  generality. 

Let  tp  G  M,  and  consider  exp(jtp)  as  a  point  on  the  complex  unit  circle.  We  define  V(tp)  as  the 
nearest  neighbor  of  exp(jtp)  in  the  set  {exp  (ju)},  v  eV,  i.e.  V(tp)  =  argmaxt^gy  cos(tP  —  v). 
Then,  our  suboptimal  quantized  beamforming  vector  p  =  Q„(h)  is  defined  by: 


\pi\  =  maxjM  :  u  <  lyA  and  u  E  U} ,  argo*  =  <  ’  ’ 

\  V(arg7/,  +  0,),  t  =  2,...,R. 

where  6q  =  V{aTgyq)  —  aigyq.  Note  that  Qu  guarantees  p  E  Cu  and  the  components  of  p  satisfy 

\Pq\  =  Ivql  =  1’  and  \yi\  -  \pi\  <  for  i  ^  q. 
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We  first  show  a  property  of  our  quantization  scheme.  Define  x  =  exp(— j6*q)p,  and  note  that 
7(x,  h)  =  7(p,h),  even  though  x  ^  Therefore,  =  7(y,  h)  -  7(p,  h)  =  7(y,  h)  - 

j(x,  h).  Thus,  for  simplicity,  we  will  use  x  for  analysis,  instead  of  p.  Similary,  we  have  \xq\  = 
\yq\  =  1,  \yi\  -  \xi\  <  for  i  ^  q.  However,  arga;,  =  argp^  -  Oq  =  mgpq  -  (V(arg|/q)  - 
argVq)  =  axgyq  as  argp^  =  V{axgyq  +  Oq)  =  V{axgyq  +  V{aigyq)  -  axgyq)  =  V(V(arg|/g))  = 
V{axgyq).  Therefore,  Xq  =  yq.  Also,  argyi  -  aigXi  =  axgyi  -  aigpi  +  6q  =  aigyi  +  6q- 
V(arg|/j  +  9q),  which  implies  |  aigyi  —  argXil  <  for  i  ^  q,hy  the  definition  of  V. 

For  the  next  step,  we  decompose  =  7(y,  h)  —  7(x,  h)  as: 

Du{h)  =  where  T>W(h)  =  7(w(*“^\h) -7(wW,h),  (27) 

i=l 

and  =  y,  =  (xi, . . .  ,Xi,  ,  pr)  for  i  =  1, . . . ,  i?  —  1,  and  =  x.  In  other 

words,  we  successively  alter  each  component  of  y,  until  we  reach  x.  First,  we  will  find  an  upper 
bound  for  D^\h),  j  G  (1, . . . ,  R},  j  ^  q.  Note  that  =  0  as  Xq  =  Pq. 

Let  6k  =  \yj\  —  \xj\  <  K~^  and  6l  =  \  aigpj  —  argXjl  <  nL~^.  Also,  pi  =  i_)_|^|2pQ  and 


A'  =  Ef=i  fiPiy^pl  =  A-  (pj  -  Xj)fjgj^j 

B'  =  1  +  Eti  =  B  -  (|%f  -  \xj\^)\9j\^pj, 


so  that  Dl^\h)  =  can  be  bounded  as 

Di  ’(h)  = - - - 

< - i - At,  '  <  2«t{.4-(s,  -  <  2\A\\i,,  -  i,||/j||ftlVft.  (28) 

Define  A*  =  Then,  we  use  the  following  inequalities  to  further  bound  (28): 


l/dv^=l/* 


i  +  \h?Po 


< -^Ai,  so  that  |A|  <  ^  ll/illt/ilVA  <  I]  (29) 


\yj  -Xj\  =  ^6\  +  2\yj\\xj\  {I  -  cos 5l)  <  \l 6\  +  \yj\\xj\6l  <  ^6\  +  5l  +  (30) 


For  (30),  the  first  inequality  can  easily  be  verified  using  the  Taylor  series  expansion  for  cos5l. 


and  we  used  |a;j|||/j|  <  1  for  the  second  inequality.  Let  us  choose  L  =  K.  Using  (29)  and  (30) 
in  (28),  iA0)(h)  <  Vl  +  7^‘^\gj\tjEE^=i  \9i\^/E^  which  is  independent  of  y  and  x.  This 
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provides  a  tractable  calculation  of  the  expectation  as  follows:  Define  Zi  =  \gi\^.  Then, 


2Vl  +  7r2 

K 

2v/rT^ 


[Zj]  Xj  + 


in  1 


R 


K 


>'!  E  E., 

^=l,  t^j 

r(2)A,  +  r(|)A|  f:  r(|)Af 

/  A 


1 1  1 


as  E^.  [Zf]  =  /q“  zfe  ^*dzi  =  r(a  +  1).  Now,  let  X'  =  max*  A*.  For  any  j  G  {1, . . . ,  R},  we  have 
Eh  where  D'^  =  2>:K-^Vi  +  1R  [l  +  (i?  -  l)T^  (|) 


which  is  independent  of  j.  Then,  using  (27)  and  the  fact  that  =  0, 

Eh  [D{h)]  <  Eh  [DM]  <  Po{R  -  l)Du- 


(31) 


Since  we  set  L  =  K,  M  =  ^  —  (K  —  1)^  .  Next,  we  show  that  —  (K  —  1)^  < 

for  any  K,  R  >  1.  We  prove  this  by  induction.  The  inequality  is  obvious  for  K  =  1.  For 
K  >2,  assume  that  it  is  true  for  k.  Then,  (k  +  1)^  —  k^  =  Y.n=i  (f)  <  Z]n=i  ^(n)  k^~"'  = 

RJ2n=o  ^  R-ik  +  1)^“\  which  shows  that  the  inequality  is  true  for  k  +  1.  Thus, 

we  have  proved  the  inequality. 

Therefore,  M  <  RK‘^^^~^\  or  equivalently,  K~^  <  2~2{fl-i)_R2{fl-i)  _  We  can  re-write  (31) 


as  Eh[T>(h)]  <  /i2  2(h-i)^  where  /i  =  2PoVl  +  ti‘^X'{R  —  l)P2(ii-i)  i  _|-  (/^  _  1)^2 
concludes  the  proof. 


.  This 
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Fig.  1:  System  block  diagram 


Fig.  2:  Quantized  power  codebook  design  for  i?  =  3  and  K  =  1,2.  The  shaded  surface  regions 
represent  the  set  of  all  feasible  power  control  vectors.  The  power  codebook  vectors  are  expressed 
by  points  on  this  region. 
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Fig.  4:  Performance  results  for  a  2-relay  network  with  unequal  individual  power  constraints: 

Po  =  Pi  =  P,  P2  =  f . 


Fig.  5:  Performance  results  for  a  3-relay  network  with  unequal  individual  power  constraints: 

Po  =  P,  Pi  =  P2  =  P3  =  f . 


Fig.  6:  BER  analysis  and  simulations  for  a  3-relay  network  with  sum-power  eonstraint 


